NASA 

TP 

1575 

c.l 


NASA Technical Paper 1575 


^LOAr-> C’ 

>.FWL 

KJRTLAND Ar i3,, : 


sN TO 

ft usy 




Determination of Noise Equivalent 
Reflectance for a Multispectral 
Scanner - A Scanner Sensitivity Study 


Daniel E. Gibbons and Richard R. Richard 


DECEMBER 1979 




TECH LIBRARY KAFB, NM 



013^7^fl 


NASA Technical Paper 1575 


Determination of Noise Equivalent 
Reflectance for a Multispectral 
Scanner - A Scanner Sensitivity Study 


Daniel E. Gibbons 
Fort Lewis College 
Durango, Colorado 

Richard R. Richard 

Lyndon B. Johnson Space Center 

Houston, Texas 


fVI/^S/\ 

National Aeronautics 
and Space Administration 

Scientific and Technical 
Information Branch 


1979 




Contents 


Section Page 

Summary 1 

Introduction 1 

Radiance 2 

Noise Equivalent Radiance 2 

Reflectance 3 

Noise Equivalent Reflectance 4 

Relationship of NER to NE Ap 4 

Determination of NE Ap 5 

Atmospheric Influence 8 

Direct Measurement of NE Ap 15 

Comparison of Measured and Calculated NE Ap 18 

Conclusion 19 

Appendix A — Meteorological Effects on Radiation 21 

Appendix B — Symbols 23 

References 26 

Tables 

Table Page 

I Solar Irradiance at Sea Level on an Area Normal to the Sun 7 

II Noise Equivalent Radiance 8 

III M^S Bandwidths 8 

IV Tabulation of NE Ap 8 

V Solar Spectral Irradiance — Proposed Standard Curve 9 

VI Standard Values of the Turbidity Coefficient ^ as a Function of Latitude and Altitude 11 

VII Summary of Available Meteorological Data for Flight Dates 13 

VIII Possible Variations in Calculated Values of NE Ap for Channel 4 of the M^S 19 


Figures 


Figure Page 

1 Schematic diagram showing coordinate system and major parameters 3 

2 Diagram of system linear in reflectance 4 

3 Ratio of diffuse irradiance to solar irradiance as a function of wavelength 12 

4 Variation of visual range with extinction coefficient at a wavelength of 0.55 micrometer . . 13 

5 Dependence of total Rayleigh optical depth on wavelength 14 

6 Aerosol optical depth as a function of visual range R^. Parameter is spectral 

wavelength X in micrometers 14 

7 Air mass m as a function of time 15 

8 Nomograph of scanner parameters for haze condition; visibility = 4 kilometers 15 

9 Nomograph of scanner parameters for clear condition; visibility = 15 kilometers 16 

10 Nomograph of scanner parameters for exceptionally clear condition; visibility = 

100 kilometers 16 

1 1 Reflectance as a function of wavelength for 5-, 8-, 16-, and 28-percent-reflectance 

panels (September 25, 1975) 16 

12 Reflectance as a function of center wavelength for 60-percent-reflectance panel 

(September 25, 1975) 16 

13 Sample printout of M^S data showing clustering method used for locating reflectance 

panel data 17 

14 Plot of scene signal voltage as a function of relative reflectance for channel 7 

(September 25, 1975) 17 

15 Plot of relative reflectance as a function of scene signal voltage for channel 3 

(September 25, 1975) 18 

16 Scanner video trace 18 

17 Geometric relationships of solar irradiance quantities used in Schwarzchild’s equation .... 21 


iv 



Determination of Noise Equivaient Ref iectance for a Muitispectrai 
Scanner — Scanner Sensitivity Study 

By Daniel E. Gibbons* and Richard R. Richard 


Summary 

Airborne remote sensor scanner data were ac- 
quired over calibrated targets on four different 
occasions, providing a measurement of scanner 
sensitivity under operational conditions. These 
measurements were then compared with predicted 
sensitivity values based on more or less standard 
calculations, one of which is performed assuming 
two sets of atmospheric conditions and a 45° Sun 
angle and two of which are based on atmospheric 
models attributable to Turner and to Angstrom. 
The comparisons show that, on the average, 
although the potential for best accuracy lies in the 
atmospherically corrected values, the assumed con- 
ditions produce good results adequate for many 
purposes. 

A nomograph was developed to provide a means 
of predicting the variations in noise equivalent re- 
flectance that are due to time of day, seasons, and 
latitude. Because of the extremes in values which 
are evident, the nomograph should prove useful to 
investigators whose requirements are stringent and 
who are faced with the use of marginally sensitive 
remote-sensing equipment. 


Introduction 

Users of airborne remote sensor scanner data are 
frequently concerned with the characterization of 
ground targets through the measurement of varia- 
tions in target reflectance. Because the variations of 
interest are often small in magnitude and are 
difficult to measure precisely, considerable interest 
is created in the definition and measurement of the 
capability of a sensor to respond to small reflec- 
tance changes. This capability, related to sensi- 
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tivity, is often described in terms of noise equiv- 
alent reflectance NE Ap, which is the minimum 
detectable variation in reflectance and is thus 
defined as the resolution of the system. 

Several methods can be used to measure or cal- 
culate this sensitivity parameter, each fraught with 
its own set of problems. Because of problems, 
mostly logistical, the usual approach is to perform a 
calculation on the basis of certain assumed target 
conditions. It is these assumed target conditions, 
which are not universally agreed upon, that are 
responsible for the variety of sensitivity numbers 
assigned to a particular instrument by various 
working groups. The possibility of improving this 
situation was one of the prime movers in this 
study, which was conducted during the summers of 
1977 and 1978. The objectives were to perform 
measurements of NE Ap and to compare these 
values with mathematically predictable values, to 
ascertain the range of values possible under normal 
operating conditions, and then to determine what 
conditions should be assumed or what approach 
should be taken to calculate the most realistic or 
meaningful values of NE Ap. 

As previously stated, data users prefer noise 
equivalent reflectance as a sensitivity parameter. 
Sensor engineers, on the other hand, prefer a 
different sensitivity parameter, noise equivalent 
radiance (NER). This preference for NER is moti- 
vated by the straightforwardness of the measure- 
ment process, which is accomplished with a min- 
imum of equipment of a variety found in most opti- 
cal laboratories. 

The first part of this report contains a limited re- 
view of terminology as an aid to those who may not 
be thoroughly familiar with, or may have forgotten, 
some of the material. Additional information on at- 
mospheric modeling is provided in appendix A. 

In compliance with the NASA’s publication pol- 
icy, the original units of measure have been con- 
verted to the equivalent value in the Systeme Inter- 



national d’Unites (SI). As an aid to the reader, the 
SI units are written first and the original units are 
written parenthetically thereafter. 


Radiance 


The detectors used in remote sensor systems 
respond to the time rate of transfer of radiant 
energy (also called radiant flux, expressed in 
watts). The most easily visualized flux measure- 
ment situation is one in which the energy-collecting 
aperture of a detector is completely filled by radia- 
tion of uniform density. Such a situation is realized 
in the use of a term called irradiance //, which is ra- 
diant flux per unit area incident on a surface given 
in watts per square centimeter. The input radiant 
power Pin (watts) is simply the product of Irra- 
diance and collecting area 


A calibrated system, then, enables the measure- 
ment of H and thereby makes the remote-sensing 
problem one of properly relating target character- 
istics to //, with due consideration of the transmit- 
ting or intervening medium. 

Target characteristics are approached through 
the concept of radiance A, which is the amount of 
power from a source of unit area radiated into a 
unit solid angle in watts per square centimeter per 
steradian. Radiance A and irradiance /fare related 
in the following way. 


H = Noj (2) 

where co is the viewing angle of the detector system 
in steradians. Thus, it can be seen, ignoring at- 
mospheric attenuation for the time being, that the 
power density H measured by the detector system 
is related to a target characteristic A through a 
detector system characteristic w. The parameters H 
and A are also defined as follows (e.g., ref. 1). 

where the quantities are defined in appendix B. 


Noise Equivalent Radiance 

The NER of a scanner system is defined as the 
amount of change in radiance required to produce a 
signal equivalent to the noise of the system. Stated 
another way, this is the quantity of radiance that 
produces a signal-to-noise ratio (S/N) of unity. The 
radiance for a signal of this .level is difficult to 
measure, and, for a linear system, a larger radiance 
change is used to determine NER. 



where A is the change in scene radiance, is the 
signal voltage proportional to the scene radiance 
change, and is the noise voltage. An input ra- 
diance change that will produce an output signal 
equal to the noise (S/N = 1) is the minimum 
detectable change in radiance level. Thus, NER is 
the minimum detectable radiance and is defined as 
the resolution of the system. 

The airborne scanners currently in use at the 
NASA Lyndon B. Johnson Space Center (JSC) use 
mechanical means to generate ground scene video 
data; that is, a rotating mirror is used to reflect 
a small element of the target scene viewed by a 
narrow-angle telescope. The image from the 
telescope is then spectrally dispersed onto detectors 
that produce voltages proportional to the collected 
spectral radiant power. The mirror rotation is con- 
tinuous so that the interior of the scanner, which 
contains visible and infrared calibration sources, is 
viewed by the detectors once for each revolution of 
the rotating scan mirror. The internal visible 
calibration source provides one means of obtaining 
the radiance change required by equation (3). In 
this application, one of the internal blackbodies is 
used to provide a zero-radiance (in the visible 
spectrum) target. Since the radiance value of the 
calibration source is known, calculation of NER for 
a linear system is a simple matter of measuring the 
change in the output signal caused by viewing first 
the blackbody and then the calibration source and, 
then, of measuring the noise signal. The NER 
calculation is completed by ratioing the signal and 
noise output voltages and dividing into the value 


2 


for radiance of the internal calibration source (from 
eq. (3)). 


AA^ 



7 ^ 2-0 

~'S/N 


A preferred way to measure NER is to use a 
more accurate laboratory source as the target. The 
reason for this preference is that, in some scanners, 
the internal calibration source optics train is not the 
same as the target-viewing optics train, and viewing 
the known radiance source through an optics train 
identical to that through which the ground target is 
viewed eliminates this possible source of error. 


Reflectance 

Reflectance is defined as the ratio of the radiant 
flux reflected by a body to that incident upon it. 
Thus, the computation or measurement of reflec- 
tance, at least in principle, is a straightforward proc- 
ess. In practice, however, the measurement prob- 
lem is complicated by such considerations as 
surface irregularities which scatter the incident 
energy in many, perhaps all, directions with 
perhaps some spatial preference. 

Other possible complications are the wavelength 
dependence of the measuring devices and the im- 
practicality of always knowing the magnitude of 
the incident radiation. These complications are ad- 
dressed and minimized as necessary, according to 
the requirements of the experimenter. However, 
generalizations or simplifying assumptions can be 
made regarding the character of targets and meas- 
urement geometry, particularly as related to reflec- 
tance and wavelength dependence. One derivation 
(ref. 2), leading to an estimation of power received 
by the detector, relates the incident power on a 
differential area dA to that reflected from the area 


z 
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FIGURE 1. — Schematic diagram showing coordinate system 
and major parameters. 

in a given direction and thereby gives rise to the 
idea of directional dependence. In general, the 
amount of radiation collected will depend on the 
direction in which the detector views the target (fig. 
1). Thus, a partial reflectance p' is defined as the 
ratio of the spectral radiance \ reflected in a given 
direction from a surface to the spectral normal irra- 
diance on an area dA. This ratio can be ex- 
pressed as 

^ ~ H 1 

where = {dP^ldA) equals the irradiance per- 
pendicular to the target and is incident power on 
dA. 

A perfectly diffusing surface is one from which 
equal amounts of energy are reflected in all direc- 
tions by an incident beam. Such a surface would ap- 
pear equally bright to a detector when viewed from 
any direction. Perfectly diffusing surfaces are as- 
sumed in this report, and they result in a simple 
relation between total reflectance and partial reflec- 
tance. 

The power reflected from the surface dP^ (fig. 1) 
can be expressed as 


dP^ = dA da (5) 


3 



where d^l is the projected solid angle defined as 


d€l = cos Q smO dd d<l> (6) 


and indicates an integration over the entire 
hemisphere above the area dA. Substituting equa- 
tions (4) and (6) into equation (5), 


7t/2 ^ 2tx 

dP^ - dA I I ^ ^ 


^ n I ft 

f / 

*^0 *'0 

( 

/ nil ^2n /dP.\ 


( 8 ) 


Because the surface is diffuse, p' is equal to a con- 
stant. Hence, equation (8) becomes 

dP^ = up’ dP^ (9) 

Finally, the ratio of dP^ /dPj is defined as the total 

reflectance p, related to the partial reflectance fac- 
tor for diffuse surfaces (eq. (9)) as given by 

P = np’ (10) 

Substituting equation (10) into equation (4), 

If the irradiance on the surface dA is at some angle 
9q measured with respect to the zenith, or a point 
directly overhead, then, 

H^I = Hj^cos6q ( 12 ) 

Finally, substituting equation (12) into equation 

( 11 ), 


N^-^H^cosd^ (13) 

If a perfectly diffusing target surface is assumed, 
equation (13) provides a useful relationship among 
target radiance, total reflectance, and irradiance. 


Noise Equivaient Reflectance 

Because a detected reflectance change signal is 
superimposed on the noise of the detecting instru- 



Reflectance 


FIGURE 2, — Diagram of system linear in reflectance. 


ment, it can be appreciated that the minimum de- 
tectable reflectance must create an output signal at 
least equal to the system noise. In fact, by defini- 
tion, the resolution of the measuring instrument is 
that amount of input which produces a signal that 
is just equal to the noise. Normally, the negligibility 
of such an input value prevents its convenient use 
and the measurement task is accomplished (with a 
linear system, fig. 2) by using a much larger input 
value and dividing by the signal-to-noise ratio. 

NE Ap = (14) 

fe) 

where Ap is the change in scene reflectance re- 
quired to produce a signal voltage of and is 
the noise voltage. Thus, NE Ap is the resolution of 
the system in terms of reflectance and is the 
minimum detectable change. 


Relationship of NER to NE Ap 

From equation (14), it can be seen that NE Ap 
depends on the change in reflectance and not on 
the absolute value of reflectance. On the other 
hand, the radiance as seen by the scanner (eq. (11)) 
will depend on both p and wavelength \. To denote 
the dependence of radiance on wavelength, the 
notation which is called spectral radiance and is 
the partial derivative of radiance with respect to 
wavelength, is used. To develop an expression for 
the relationship between spectral radiance and the 
combination of p and X, a linear approximation 
may be used for A^;^(p,\) about some values PqAo^ 


j 
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where the change in spectral radiance is denoted as 
AA/,. 

dN. dN, 

where 

/w, = - N^{Po,\) 

Dividing each side of equation (15) by the signal- 
to-noise ratio produces the spectral noise equiv- 
alent radiance (NER;^). 


Integrating equation (20) with respect to X results 
in 


ner = n^ 



cos 0Q C^X + 



S^dX 


( 21 ) 


Then, if the factor in brackets in equation (21) is 
denoted as , the total target irradiance, one has 

NER = -///NE Ap (22) 

7T ' 


S/N 


NER, 


Equation (22) is a very useful relationship in deter- 
mining NE Ap. 


dN. Ap dN. AX 

S/N S/N (1^1 

The scene or target spectral radiance can be re- 
lated to the ground irradiance (ref. 3) by 

+ (17) 


Determination of NE Ap 

Several means of determining NE Ap for a scan- 
ner system are available; three methods are dis- 
cussed in the following paragraphs. The first of 
these is the most direct approach to measurement; 
it is based on the assumption of a linear system and 
involves the scanning of two targets of known 
reflectance. 


where is the diffuse or sky irradiance and is 
the path radiance due to scattering of radiation into 
the beam from surrounding ground targets. Taking 
the partial derivatives of equation (17) with respect 
to p and X results in 






dp 




cos 0Q + 


dN. 


t\ 


dX 


71 9X 


dN 


r\ 


ax 


(18) 

(19) 


As a first approximation, it is assumed that 



dH, 




ax 


0 ; 


dX 


= 0 


NEAp = 


P2 Pi 



(23) 


where and V 2 are the output voltages produced 
by the scanner when aimed at targets 1 and 2, and 
is the noise voltage of the system. Note that 
equation (23) is equivalent to equation (14). 

A second means of calculating NE Ap can be 
derived from equation (13). With the realization 
that the two variables are the radiance and the 
reflectance p, write 

AA^, = —H. cos0„ 

A. Tj- A. U 

where the deltas represent changes in the quan- 
tities. Dividing by S/N, 


These approximations are defended on the basis 
that the variation in the given quantities is small 
with respect to wavelength. Then, substituting 
equations (18) and (19) into equation (16), noting 
that Ap/(S/N) = NE Ap, 


S/N 


NER^ 


Ap 

S/N ’ w 


( 20 ) 


= NEAp 


cos 0Q 

7T 
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Rearranging, 


NEAp = 


7T*NER^ 
cos 0Q 


(24) 


It can be seen that using equation (24) requires that 
the solar irradiance, must be known through 
measurement or calculation or through assump- 
tions involving longitude, atmospheric scattering, 
transmission, time of day, etc. 

The third method of determining NE Ap is a 
real-time measurement that uses the internal cali- 
bration source of the scanner and a fictitious (but 
hopefully realistic) solar irradiance value. The 
method involves a computation of the equivalent 
reflectance of the internal calibration source.^ 



2.00, an atmospheric pressure of 101.3 kN/m^ (760 
torr), and a water-vapor pressure of 2.7 kN/m^ (20 
torr) are assumed with no provisions for sky or 
diffuse irradiance. These assumptions bring the 
measured and calculated numbers into fair agree- 
ment. Combining equations (25) and (26), one has 


NE Ap = 


7riVj^(cal) 
0.707^^ * 


1 

S/N 


TfNER^ 

0.707//^ 


(27) 


As calculated from equation (27), NE Ap is 
designated as the P. Moon value, the one most used 
at JSC for the determination of NE Ap for any par- 
ticular scanner. To calculate NE Ap from equation 
(27), a numerical integration must be performed 
and this removes the subscript from the two terms. 


where Pg^ is equivalent scene reflectance. In this 
expression, the energy from the calibration lamp is 
treated as though it were an input from a ground 
target caused by an external energy source, such as 
the Sun. To have physical significance, Pg^ ob- 
viously must represent a ratio (of energy reflected 
to energy incident) after the manner of equation 
(4). Substituting equations (10) and (12) into equa- 
tion (4), 


7rTV^(cal) 

//j^COS0Q 


(26) 


xTtr A - ^*NER 

^ ^ AX 

tt-NER 

0.707//p 


(28) 


To illustrate the use of equation (28), NE Ap 
will be calculated for the data of September 25, 
1975, for channel 2 of the modular multispectral 
scanner (M^S), the scanner used to acquire the data 
for this project. The following table will make the 
procedure more easily understood. 


where AJ^(cal) is the calibration source spectral ra- 
diance and is an assumed target irradiance pro- 
duced by incident solar radiation; is the irra- 
diance that would produce a signal Ag^| after reflec- 
tion from a fictitious target of reflectivity Pgq, 
Thus, the usefulness of NE Ap as a sensitivity 
parameter as determined by equation (25) with the 
use of equation (26) is strongly dependent on the 
degree to which approximates a true or realistic 
operational condition. In the JSC Engineering and 
Development Directorate, a nadir angle of 45° is 
assumed and is taken from a table compiled by 
P. Moon (table I) in which a relative air mass of 


*F. Fulton, Data Processing Requirements for Engineering 
Evaluation of Bendix MSDS, Lockheed Electronics Company, 
July 1971. 


Item 

Value 

Source 

NER 

4.4X10”' W.cm“^sr' 

Table II 

Center wavelength 

0.465 micrometer 

Table III 

Bandwidth 

0.03 micrometer 

Table III 

Solar irradiance 

0.1110 

Table I 


(P. Moon) 


Then, 

ttX 4.4X10“^ 

0.707X 0.1110X 0.03 

= 5.86 X 10“3 

and, expressed as a percentage, NE Ap for channel 
2 of the M^S on September 25, 1975, was 0.59 per- 
cent. (Other NE Ap values are listed in table IV.) 
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Table I. — Solar Irradiance at Sea Level on an Area Normal to the Sun' 


[m=2;//o= 1322W.m"’“] 




X , 

Hax^ 

X , 

Hax- 

X , fim 

Hax- 






W*m ‘ 


lV>m 

0.301 

0.177 

0.68 

1149 

1.23 

420 

1.78 

18.4 

.302 

.342 

.69 

978 

1.24 

387 

1.79 

5.70 

.303 

.647 

.70 

1108 

1.25 

328 

1.80 

.920 

.304 

1.16 

.71 

1070 

1.26 

311 

1.81 

— 

.305 

1.91 

.72 

832 

1.27 

381 

1.82 

— 

.306 

2.89 

.73 

965 

1.28 

382 

1.83 



.307 

4.15 

.74 

1041 

1.29 

346 

1.84 

— 

.308 

6.11 

.75 

867 

1.30 

264 

1.85 

— 

.309 

8.38 

.76 

566 

1.31 

208 

1.86 

— 

.310 

11.0 

.77 

968 

1.32 

168 

1.87 

— 

.311 

13.9 

.78 

907 

1.33 

115 

1.88 



.312 

17.2 

.79 

923 

1.34 

58.1 

1.89 

— 

.313 

21.0 

.80 

857 

1.35 

18.1 

1.90 

— 

.314 

25.4 

.81 

698 

1.36 

.660 

1.91 

.705 

.315 

30.0 

.82 

801 

1.37 

— 

1.92 

2.34 

.316 

34.8 

.83 

863 

1.38 



1.93 

3.68 

.317 

39.8 

.84 

858 

1.39 

— 

1.94 

5.30 

.318 

44.9 

.85 

839 

1.40 

— 

1.95 

17.7 

.319 

49.5 

.86 

813 

1.41 

1.91 

1.96 

31.7 

.32 

54.0 

.87 

798 

1.42 

3,72 

1.97 

37.7 

.33 

101 

.88 

614 

1.43 

7,53 

1.98 

22.6 

.34 

151 

.89 

517 

1.44 

13,7 

1.99 

1.58 

.35 

188 

.90 

480 

1.45 

23,8 

2.00 

2.66 

.36 

233 

.91 

375 

1.46 

30.5 

2.01 

19.5 

.37 

279 

.92 

258 

1.47 

45.1 

2.02 

47.6 

.38 

336 

.93 

169 

1.48 

83.7 

2.03 

55.4 

.39 

397 

.94 

278 

1.49 

128 

2.04 

54.7 

.40 

470 

.95 

487 

1.50 

157 

2.05 

38.3 

.41 

672 

.96 

584 

1.51 

187 

2.06 

56-2 

.42 

733 

.97 

633 

1.52 

209 

2.07 

77.0 

.43 

787 

.98 

645 

1.53 

217 

2.08 

88.0 

.44 

911 

.99 

643 

1.54 

226 

2.09 

86.8 

.45 

1006 

1.00 

630 

1.55 

221 

2.10 

85.6 

.46 

1080 

1.01 

620 

1.56 

217 

2.11 

84.4 

.47 

1138 

1.02 

610 

1.57 

213 

2.12 

83.2 

.48 

1183 

1.03 

601 

1.58 

209 

2.13 

20.7 

.49 

1210 

1.04 

592 

1.59 

205 

2.14 

— 

.50 

1215 

1.05 

551 

1.60 

202 



.51 

1206 

1.06 

526 

1.61 

198 



.52 

1199 

1.07 

519 

1.62 

194 



.53 

1188 

1.08 

512 

1.63 

189 



.54 

1198 

1.09 

514 

1.64 

184 



.55 

1190 

1.10 

252 

1.65 

173 



.56 

1182 

1.11 

126 

1.66 

163 



.57 

1178 

1.12 

69.9 

1.67 

159 



.58 

1168 

U 3 

98.3 

1.68 

145 



.59 

1161 

1.14 

164 

1.69 

139 



.60 

1167 

1.15 

216 

1.70 

132 



.61 

1168 

1.16 

271 

1.71 

124 



.62 

1165 

1.17 

328 

1.72 

115 



.63 

1176 

1.18 

346 

1.73 

105 



.64 

1175 

1.19 

344 

1.74 

97.1 



.65 

1173 

1.20 

373 

1.75 

80.2 



.66 

1166 

1.21 

402 

1.76 

58.9 



.67 

1160 

1-22 

431 

1.77 

38.8 




®From reference 4. See appendix B for definition of symbois. 
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Table II. — Noise Equivalent Radiance 


Table IV . — Tabulation of NE Ap 


[Watts per square centimeter per steradian] 


[Percent] 


Channel 

Group I 

Group II 

Group III 

Group IV 

1 

_ 

22.4 X 10 ”* 

12.9 X 10 "* 

13.7 X 10 "‘ 

2 

4.6 X 10 "° 

7.1 

— 

4.7 

3 

4.1 

5.01 

4.36 

3.70 

4 

4.1 

4.1 

3.55 

4.20 

5 

— 

3.3 

3.7 

3.9 

6 

4.0 

3.6 

3.6 

4.1 

7 

3.7 

4.3 

3.6 

3.8 

8 

3.4 

6.7 

3.4 

3.0 

9 

4.1 

11.0 

3.3 

1.1 

10 

4.4 

7.7 

8.09 

11.43 


Table III. — M^S Bandwidths 
[Micrometers] 

Channel 

Bandpass 

Bandwidth 

Center 




wavelength 

1 

0.41 to 0.44 

0.03 

0.425 

2 

0.4510 0.48 

.03 

.465 

3 

0.49 to 0.53 

.04 

.51 

4 

0.53 to 0.57 

.04 

.55 

5 

0.57 to 0.61 

.04 

.59 

6 

0.61 to 0.65 

.04 

.63 

7 

0.65 to 0.69 

.04 

.67 

8 

0.69 to 0.74 

.05 

.71 

9 

0.76 to 0.86 

.10 

.81 

10 

0.95 to 1.03 

.08 

.99 

11 

7.89 to 11.64 

3.75 

9.77 


Atmospheric Influence 

Equation (22) shows the relationship of NE Ap 
to the scene irradiance and the scanner charac- 
teristic, NER. To make use of this relationship, one 
must be able to predict or calculate with some 
degree of accuracy, and it is here one must consider 
the effect of the atmosphere. Several atmospheric 
models are available; one is attributable to 
Angstrom (ref. 5) and has the following form. 




m (a^+a) 

[ R al w o 


(29) 


where 9q is the zenith Sun angle, is the extrater- 
restrial spectral irradiance (from Thekaekara, table 
V), and m is air mass and is equal to secant 9 for 
angles between 0° and 60°. The other terms are 


Channel 

Group I Group II 

Group III 

Group IV 

Measured value 

1 



1.65 

0.74 

0.68 

2 

0.50 

.40 

— 

.19 

3 

.41 

.34 

.16 

.16 

4 

.49 

.27 

.27 

.21 

5 

— 

.27 

.27 

.23 

6 

.51 

.32 

.25 

.23 

7 

.53 

.33 

.25 

.21 

8 

.55 

.36 

.24 

.18 

9 

.39 

.37 

.16 

.36 

10 

.66 

.55 

2.30 

3.25 



P. Moon value 



1 



4.37 

2.51 

2.67 

2 

0.59 

.95 

— 

.63 

3 

.38 

.46 

.40 

.34 

4 

.38 

.38 

.33 

.39 

5 

— 

.32 

.35 

.37 

6 

.38 

.34 

.34 

.39 

7 

.35 

.41 

.34 

.36 

8 

.28 

.55 

.28 

.25 

9 

.26 

.70 

.21 

.49 

10 

.38 

.67 

.70 

.99 

o 

Angstrom value 

1 



1.57 

0.94 

1.00 

2 

0.31 

.39 

— 

.26 

3 

.20 

.21 

.19 

.16 

4 

.23 

.18 

.16 

.19 

5 

— 

.15 

.17 

.18 

6 

.24 

.17 

.18 

.20 

7 

.25 

.23 

.19 

.21 

8 

.20 

.31 

.16 

.15 

9 

.16 

.34 

.11 

.25 

10 

.34 

.45 

.51 

.71 

Nomograph value 

1 



3.88 

0.80 

0.85 

2 

0.62 

.87 

— 

.24 

3 

.38 

.44 

.17 

.14 

4 

.39 

.35 

.15 

.18 

5 

— 

.27 

.15 

.16 

6 

.38 

.31 

.16 

.18 

7 

.36 

.38 

.17 

.18 

8 

.29 

.50 

.14 

.13 

9 

.21 

.51 

.09 

.20 

10 

.41 

.60 

.38 

.53 
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Table V. — Solar Spectral Irradiance — Proposed Standard Curve^ 
[Solar constant = 1353 W»m"^l 


X, fxm 

Efx;, 

W'.m ' 


D(O-K), 

percent 

X. fxm 

E^X;. 

E(O-K). 

percent 

0.120 

0.100 

0.0059992 

0.00044 

0.375 

1157 

89.0641 

6.582 

.140 

.030 

.0072993 

.00053 

.380 

1120 

94.7566 

7.003 

.150 

.07 

.00780 

.00057 

.385 

1098 

100.3016 

7.413 

.160 

.23 

.00930 

.00068 

.390 

1098 

105.7916 

7.819 

.170 

.63 

.01360 

.00100 

.395 

1189 

111.5091 

8.241 

.180 

1.25 

.02300 

.00169 

.400 

1429 

118.0541 

8.725 

.190 

2.71 

.04280 

.00316 

.405 

1644 

125.7366 

9.293 

.200 

10.7 

.10985 

.00811 

.410 

1751 

134.2241 

9.920 

.210 

22.9 

.27785 

.02053 

.415 

1774 

143.0366 

10.571 

.220 

57.5 

.67985 

.05024 

.420 

1747 

151 8391 

11 222 

.225 

64.9 

.98585 

.0728 

.425 

1693 

160.4391 

11 858 

.230 

66.7 

1.31485 

.0971 

.430 

1539 

168 7691 

12.473 

.235 

59.3 

1.62985 

.1204 

.435 

1563 

177 0241 

13 083 

.240 

63.0 

1.93560 

.1430 

.440 

1810 

185.7066 

13 725 

.245 

72.3 

2.27385 

.1680 

.445 

1922 

195 0356 

14.415 

.250 

70.4 

2.63060 

.1944 

.450 

2006 

204 8566 

15 140 

.255 

104 

3.06660 

.2266 

.455 

2057 

215 0141 

15.891 

.260 

130 

3.65160 

.2698 

.460 

2066 

225 3216 

16 653 

.265 

185 

4.43910 

.3280 

.465 

2049 

235 6066 

17413 

.270 

232 

5.48160 

.4051 

.470 

2033 

245 8091 

18 167 

.275 

204 

6.5716 

.4857 

.475 

2044 

256 001 

18 921 

.280 

222 

7.6366 

.5644 

.480 

2074 

266 296 

19 681 

.285 

315 

8.9791 

.6636 

.485 

1976 

276.421 

20.430 

.290 

482 

10.9716 

.8109 

.490 

1950 

286 236 

21.155 

.295 

584 

13.6366 

1.0078 

.495 

1960 

296.011 

21.878 

.300 

514 

16.3816 

1.2107 

.500 

1942 

305 766 

22 599 

.305 

603 

19.1741 

1.4171 

.505 

1923 

315 421 

23 312 

.310 

689 

22.4041 

1.6558 

.510 

1982 

324 926 

24 015 

.315 

764 

26.0366 

1.9243 

.515 

1933 

334 214 

24 701 

.320 

830 

30.0216 

2.2188 

.520 

1933 

343 379 

25 379 

.325 

975 

34.5341 

2.552 

.525 

1852 

352 591 

26 059 

.330 

1059 

39.6191 

2.928 

.530 

1842 

361 826 

26 742 

.335 

1081 

44.9691 

3.323 

.535 

1819 

370 976 

27 418 

.340 

1074 

50.3566 

3.721 

.540 

1783 

379 979 

28 084 

.345 

1069 

55.7141 

4.117 

.545 

1754 

388 821 

28 737 

.350 

1093 

61.1191 

4.517 

.550 

1725 

397 519 

29 380 

.355 

1083 

66.5591 

4.919 

.555 

1720 

406.131 

30 017 

.360 

1068 

71.9366 

5.316 

.560 

1695 

414 669 

30 648 

.365 

1132 

77.4366 

5.723 

.565 

1705 

423.169 

31.276 

.370 

1181 

83.2191 

6.150 

.570 

1712 

431 711 

31 907 


^From reference 6 See appendix B for definiiion of symbols 


given in appendix B but some discussion is in 
order. The quantity is the Rayleigh extinction 
coefficient (or scattering coefficient) and it pro- 
vides a measure of the reflective interaction be- 
tween the solar radiant energy and the molecules of 


the gaseous atmosphere. This coefficient is wave- 
length dependent and is given by 

a^(X) = 0.00889X“'^-°^ (30) 
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Table V . — Continued 


X, fim 

E(K). 

E(O-k), 

D(O-k). 

percent 

X, /xm 

EjXX 

W*m 

E(O-k). 

1 IJ/ 

D(()-K), 

percent 

0.575 

1719 

440.289 

32.541 

0.94 

847 

892.08 

65.934 

.580 

1715 

448.874 

33.176 

.95 

837 

900.50 

66.556 

.585 

1712 

457.441 

33.809 

.96 

820 

908.79 

67.168 

.590 

1700 

465.971 

34.439 

.97 

803 

916.90 

67.768 

.595 

1682 

474.426 

35.064 

.98 

785 

924.84 

68.355 

.600 

1666 

482.796 

35.683 

.99 

767 

932.60 

68.928 

.605 

1647 

491.079 

36.295 

1.00 

748 

940.18 

69.488 

.610 

1635 

499.284 

36.902 

1.05 

668 

975.58 

72.105 

.620 

1602 

515.469 

38.098 

1.10 

593 

1007.10 

74.435 

.630 

1570 

531.329 

39.270 

1.15 

535 

1035.30 

76.519 

.64 

1544 

546.899 

40.421 

1.20 

485 

1060.80 

78.404 

.65 

1511 

562.174 

41.550 

1.25 

438 

1083.88 

80.101 

.66 

1486 

577.159 

42.557 

1.30 

397 

1104.75 

81.652 

.67 

1456 

591.869 

43.744 

1.35 

358 

1123.63 

83.047 

.68 

1427 

606.284 

44.810 

1.40 

337 

1141.00 

84.331 

.69 

1402 

620.429 

45.855 

1.45 

312 

1157.23 

85.530 

.70 

1369 

634.284 

46.879 

1.50 

288 

1172.23 

86.639 

.71 

1344 

647.849 

47.882 

1.55 

267 

1186.10 

87.665 

.72 

1314 

661.139 

48.864 

1.60 

245 

1198.40 

88.611 

.73 

1290 

674.159 

49.826 

1.65 

223 

1210.60 

89.475 

.74 

1260 

686.909 

50.769 

1.70 

202 

1221.23 

90.261 

.75 

1235 

699.384 

51.691 

1.75 

180 

1230.78 

90.967 

.76 

1211 

711.614 

52.595 

1.80 

159 

1239.25 

91.593 

.77 

1185 

723.594 

53.480 

1.85 

142 

1246.78 

92.149 

.78 

1159 

735.314 

54.346 

1.90 

126 

1253.48 

92.644 

.79 

1134 

746.779 

55.194 

1.95 

114 

1259.48 

93.088 

.80 

1109 

757.994 

56.023 

2.00 

103 

1264.90 

93.489 

.81 

1085 

768.966 

56.834 

2.10 

90 

1274.55 

94.202 

.82 

1060 

779.694 

57.627 

2.20 

79 

1283.00 

94.826 

.83 

1036 

790.174 

58.401 

2.30 

69 

1290.40 

95.373 

.84 

1013 

800.419 

59.158 

2.4 

62.0 

1296.95 

95.8580 

.85 

990 

810.434 

59.899 

2.5 

55.0 

1302.80 

96.2903 

.86 

968 

820.224 

60.622 

2.6 

48.0 

1307.95 

96.6710 

.87 

947 

829.793 

61.330 

2.7 

43.0 

1312.50 

97.0073 

.88 

926 

839.164 

62.022 

2.8 

39.0 

1316.60 

97.3103 

.89 

908 

848.334 

62.700 

2.9 

35.0 

1320.30 

97,5838 

.90 

891 

857.329 

63.365 

3.0 

31.0 

1323.60 

97.8277 

.91 

880 

866.184 

64.019 

3.1 

26.0 

1326.45 

98.0383 

.92 

869 

874.929 

64.665 

3.2 

22.6 

1328.88 

98.2179 

.93 

858 

883.564 

65.304 

3.3 

19.2 

1330.97 

98.3724 


where X is given in micrometers. The quantity is 
an aerosol attenuation coefficient due to dust, haze, 
etc., and is approximated by 

fl^(X) = 2.303£(2X)"“ (31) 

which introduces two more quantities, B and a. 
The term 5 is called a turbidity constant, the values 
for which range from 0.05 to 0.2 (table VI); a takes 


on a value 0 < a ^ 4, depending on particle size, 
and is given the value 1.0 for light haze and the 
value 1.5 for clear conditions. The quantity n„,is an 
attenuation coefficient that varies with latitude and 
wavelength but, for a water-vapor pressure of 2.7 
kN/m^ (20 torr), is 0.011 for 0.4 < X ^ 0.7 mi- 
crometer and 0.18 for 0.6 ^ X ^ 1.0 micrometer. 
The quantity S^, diffuse sky irradiance, is highly 
variable but can be approximated from experimen- 
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Table V, — Concluded 


X,^m E(X;, _ \^( 0 -\), 

W» m~^*fJLmT ^ IV^ m~^ percent 


Table VI. — Standard Values of the Turbidity 
Coefficient % as a Function of Latitude and Altitude 


Altitude 

Level 

Latitude 

(air 

pressure). 


70^ N 60° N 45° N 30° N 0‘ 

kN/m^ 

(mbar) 




100(1000) Angstrom 0.050 0.061 


3.4 

16.6 

3.5 

14.6 

3.6 

13.5 

3.7 

12.3 

3.8 

11.1 

3.9 

10.3 

4.0 

9.5 

4.1 

8.7 

4.2 

7.8 

4.3 

7.1 

4.4 

6.50 

4.5 

5.90 

4.6 

5.30 

4.7 

4.80 

4.8 

4.50 

4.9 

4.10 

5.0 

3.83 

6.0 

1.75 

7.0 

.99 

8.0 

.60 

9.0 

.380 

10.0 

.250 

il.O 

.170 

12.0 

.120 

13.0 

.087 

14.0 

.055 

15.0 

.049 

16.0 

.038 

17.0 

.031 

18.0 

.024 

19.0 

.02000 

20.0 

.01600 

25.0 

.00610 

30.0 

.00300 

35.0 

.00160 

40.0 

.00094 

50.0 

.00032 

60.0 

.00019 

80.0 

.00007 

100.0 

.00001 


1332.76 

98.5047 

1334.32 

98.6200 

1335.73 

98.7238 

1337.02 

98.8192 

1338.19 

98.9056 

1339.26 

98.9847 

1340.25 

99.0579 

1341.16 

99.1252 

1341.98 

99.1861 

1342.73 

99.2412 

1343.4141 

99.291507 

1344.0341 

99.337331 

1344.5941 

99.378721 

1345.0991 

99.416045 

1345.5641 

99.450413 

1345.9941 

99.482195 

1346.3906 

99.511500 

1349.1806 

99.717708 

1350.5506 

99.818965 

1351.3456 

99.877723 

1351.8356 

99.913939 

1352.1506 

99.937221 

1352.3606 

99.952742 

1352.5056 

99.963459 

1352.6091 

99.971108 

1352.6801 

99.976356 

1352.7321 

99.980199 

1352.7756 

99.983414 

1352.8101 

99.985964 

1352.8376 

99.987997 

1352.8526 

99.985623 

1352.8776 

99.990953 

1352.9228 

99.995037 

1352.9556 

99.996718 

1352.9671 

99.997568 

1352.9734 

99.998037 

1352.9800 

99.998525 

1352.9829 

99.998736 

1352.9855 

99.998928 

1352.9865 

99.999002 

1353.0000 

100.000000 




Lower limit 

.010 

.010 



Upper limit 

.100 

.100 

90 

(900) 

Angstrom 

.030 

.037 



Lower limit 

.010 

.010 



Upper limit 

.050 

.100 

80 

(800) 

Angstrom 

.018 

.022 



Lower limit 

.000 

.000 



Upper limit 

.050 

.050 


1.082 

0.104 

0.125 

0.104 

0.082 

0.061 

.050 

.050 

.050 

.050 

.050 

.010 

.200 

.400 

.400 

.400 

.200 

.100 

.050 

.063 

.076 

.063 

.050 

.037 

.010 

.010 

.010 

.010 

.010 

.010 

.100 

.200 

.400 

.200 

.100 

.100 

.030 

.037 

.046 

.037 

.030 

.022 

.010 

.010 

.010 

.010 

.010 

.000 

.100 

.200 

.200 

.200 

.100 

.050 

.018 

.023 

.028 

.023 

.018 

.014 

.000 

.010 

.010 

.010 

.000 

.000 

.050 

.100 

.200 

.100 

.050 

.050 


1 000.0 .00000 


tal evidence which shows that SJH>^ = 0.35 for 
wavelengths between 0.4 and 0.6 micrometer, and 
SJH^ = 0.1 for wavelengths between 0.6 and 1.0 
micrometer (ref. 5). 

The wavelength dependency of the exponents in 
equation (29) results in a formidable expression. 
However, simplifications are possible in a measure- 


70 (700) Angstrom .011 .014 

Lower limit .000 .000 

Upper limit .020 .050 


ment situation because the relatively narrow bands 
used in a scanner system permit the use of approx- 
imations within a band. Thus, to find the scene irra- 
diance in a band corresponding to the detector 
responses, an integration of equation (29) would be 
performed in the following way. 


/ 


H^d\ = H, 




cos 0 










(32) 


where 




-L 




and 


= 


/ 


S^d\ 


and the (X^)’s are the center wavelengths of the 
band passes. Finally, NE Ap is computed for a 
given scanner and considering the atmospheric 
attenuation. 


NEAp = 


7t«NER 

H. 


(33) 
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where is the value taken from equation (32). 
Another atmospheric model, which deals with the 
terms optical depth and visibility, is discussed in 
the following paragraphs and in appendix A. 

To illustrate the use of equations (32) and (33), 
NE Ap will be computed using M^S data acquired 
on May 15, 1978. The method is best illustrated 
using the following table. 


(channel 4) then has a value of 0,33 for 
Then, 


= H^X 0.33 


= 1725 X X 0.33 

= 5.69 X 10“2 


and 


Item 

Value 

Source 

5, = J's^dX^S^A-K 

e 

26.06“ 

Table VII 

cos B 

0.90 

Computation 

= 5.69 X 10^2 AX 

m 

1.11 

m = secant 9 


0.55 micrometer 

Table III (channel 4) 



0.10 

= 0.00889X”"“ 

= 2.28 X 10“^ 

B 

0.082 

Table VI, latitude 45' 


OL 

1.5 

Assumption 

Then, from equation (32), 


0.16 

2.3035(2X)-“ 


K 

0.011 

Text 

1725 

AX 

0.04 micrometer 

Table III (channel 4) 



1725xl0“‘'w.cm~^/xm“‘ 

Table V (0.55 micrometer) 

+ 2.28 X 10^^ 


2.28xlO"^W.cm“'‘ 

Text 


NER 

3.55xl0"'’w.cm“^sr“' 

Table II 

= 6.88 X 10~^ 


According to Robinson (ref. 5), 5 is approxi- 
mated from experimental evidence which shows 
that S^/H^ ^ 0.35 for wavelengths between 0,4 and 
0,6 micrometer and ^0.1 for wavelengths 

between 0.6 and 1.0 micrometer. A curve is drawn 
(fig. 3) to approximate this condition and remove 
the discontinuity. The 0.55-micrometer band 


Then, substituting into equation (33), 
7t*NER 


NEAp = 




7T X 3.55 X 10“^ 


6.88 X 10 


-3 



FIGURE 3. — ^Ratio of diffuse irradiance to solar irradiance as 
a function of wavelength. 


= 0.16 percent 

Other values of Angstrom’s NE Ap are given in ta- 
ble IV. 

Assuming that the Angstrom or any other 
model describes at least the average properties of 
the atmosphere, it is possible to develop a range of 
values for NE Ap under given atmospheric condi- 
tions; and atmospheric conditions are best de- 
scribed, in a practical sense, in terms of visibility or 
visual range. Visual range is related to aerosol or 
haze scattering and is a quantity quoted by the Na- 
tional Weather Service as a measure of the gross 
features of atmospheric conditions. These have 
been classified and correspond to the attenuation of 
the atmosphere at 0.55 micrometer as haze, clear. 
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Table VII. — Summary of A vailabie Meteorological Data for Flight Dates 


Item 


Data groups— flight dates 



Group I — 

Group II — 

Group III — 

Group IV — 


Sept. 25, 1975 

July 8, 1976 

May 15, 1978 

June 2, 1978 

Latitude, deg N 

44.563 

44.563 

44.563 

44.563 

Longitude, deg W 

98.908 

98.908 

98.908 

98.908 

Altitude, m (ft) 

451 (1480) 

451 (1480) 

451 (1480) 

451 (1480) 

Local time 

1346 

1343 

1337 

1228 

Sun angle, deg 

44.50 

22.50 

26.06 

26.96 

Solar declination, deg 

-0.520139 

22.4862 

18.7265 

22.1083 

Relative air mass 

1.402 

1.082 

1.079 

1.122 

Cloud cover, percent 

NA^ 

5 

20 

0 

Relative humidity, percent 

NA 

30 

54 

35 

Barometric pressure, kN/m^ (in. Hg) 

NA 

96.17 (28.40) 

96.24 (28.42) 

96.85 (28.60) 

General conditions 

NA 

Haze, dust 

Clear 

Clear 


®Not available. 


etc. The optical depth of the atmosphere to haze at- 
tenuation is then related to visual range, which is 
the distance at which the ratio of the apparent 
contrast of an object viewed with background radia- 
tion, viewed by an observer at range R, to the in- 
herent contrast at zero range is equal to 2 percent. 
Figure 4 shows the relationship of the terms 
“clear,” “very clear,” “haze,” etc., with daylight 
visibility range in kilometers. Atmospheric models 
that are used to correct remote-sensing data have 
been developed which make use of the visual range 
quantity. One such model, attributable to Turner 
(ref. 7), has the form 


- m { T n + T \ 

= H^cose^e ^ ^ (34) 

There is a great deal of similarity in the ap- 
pearance of equation (34), Turner’s model, and 
equation (29), Angstrom’s model. Turner’s model 
does not include a water absorption quantity as 
does Angstrom’s model as in the term a^. Since it is 
primarily the visible region of the spectrum that is 
involved, there is no inconsistency because a^ is 
very small below 0.72 micrometer, and Angstrom’s 
model is simply more general. The attractiveness of 
equation (34) lies in the use which is made of the 
actual condition, visual range. The exponent in 
equation (34) introduces the terms and r^, which 
are the total Rayleigh optical depth (fig. 5) and the 
aerosol optical depth (fig. 6), respectively. The 
subscript t is used on the sky irradiance term to 



FIGURE 4. — Variation of visual range with extinction coeffi- 
cient at a wavelength of 0.55 micrometer. 


make a distinction from the term in equation 
(34); S^is found from the relationship (ref. 7) 


k 

m 


(35) 
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Wavelength, |im 

FIGURE 5. — ^Dependence of total Rayleigh optical depth on 

wavelength. 

and for the data in this report, 



Visual range, km 


FIGURE 6. — Aerosol optical depth as a function of visual 
range Parameter is spectral wavelength X in micrometers. 


k = qm - 0.68 


where 




q = e 


(36) 


(37) 


The value 0.68 in equation (36) is an empirically 
derived constant related primarily to the scanner 
characteristics. From equation (33), 


NE Ap = ttNER (33) 

Using equations (34) and (35), substituting in equa- 
tion (33), and rearranging terms. 


NE Ap 


mq'» 

_NER _ 

_ 

\ ^b 


(38) 


An inspection reveals that the left side of equation 
(38) consists of scanner terms and the right side 
relates to physical conditions, which in turn relate 
to the latitude of the observation, the time of day, 
and the state of the atmosphere, etc. With equa- 
tions (35) and (36), equation (38) becomes 


NE Ap 


.NER. 



m 


1 - 


mq 

■^8 

m 


(39) 


Equation (39) forms the basis of nomographs 
wherein one may resolve scanner system sen- 
sitivity on the basis of the conditions relating to the 
data situation. The value for m may be calculated or 
may be taken from the graph in figure 7, which is a 
plot of air mass as a function of time T in hours 
before or after noon for two latitudes and three 
seasons of the year. This value of m is used with 
figures 8, 9, and 10 to find a value, a unitless num- 
ber, for the left side of equation (39). Then, know- 
ing the scanner NER and the extraterrestrial irra- 
diance from table V for the center wavelength of a 
particular scanner band, one may calculate NE Ap. 
Following is a step-by-step procedure for calculat- 
ing NE Ap. 

1. Using figure 7, estimate the air mass for a 
given latitude and time of day. The air mass may 
also be calculated as indicated in appendix A since 
m = 1/cos 9. 

2. Using figure 4 and an estimate of the visi- 
bility conditions, find the daylight visibility range. 

3. Using the appropriate figure (8, 9, or 10, cor- 
relative to haze, clear, or exceptionally clear) for 
the ratio of system parameters, estimate the NE Ap 
value for that channel given tt, //^, and NER. 

As an example, consider the data obtained with 
the M^S for channel 2 on June 2, 1978: (1) = 

0.47 micrometer (table III), (2) AX = 0.03 micro- 
meter (table III), (3) NER = 4.7 x 10"® 
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W.cm^2.sj.-i (table II), (4) A\ = 2033 

X 10"^ X 0,03 = 6,1 X 10-3 (tables III and V), (5) 
time of day = 12:30 p,m. = 0.5 hour after noon (ta- 
ble VII), (6) latitude = 44.5° N (table VII), and (7) 
air mass = 1.122 (table VII). From figure 7 for 0.5 
hour after noon, an air mass of approximately 1.0 is 
found. (Summer solstice occurs June 21 or 22.) 
First, considering a hazy day, a value of 1.1 is ob- 
tained for the scanner constants (fig, 8). Thus, 


NE Ap = 


tt-NER X 


H 


1.1 


7T X 4.7 X 10~^ X 1.1 
6.1 X 10“3 


= 0.27 percent 

The figure 0.27 percent, calculated assuming haze, 
is the worst condition. Assuming clear conditions 
(visibility equals 15 kilometers), the best value of 
0.23 percent is obtained. This compares with the 
experimental or measured value of 0.19 percent. 
Hence, the expected value at this location and time 
would be about 0.25 percent regardless of the at- 
mospheric state. This constancy is due, of course, 
to the optimum time of the particular data run. 



FIGURE 8. — Nomograph of scanner parameters for haze con- 
dition; visibility = 4 kilometers. 


Direct Measurement of NE Ap 

One of the major activities of the JSC Earth 
Resources Program is the Large Area Crop Inven- 
tory Experiment (LACIE). The LACIE is a long- 
term project designed to develop remote-sensing 
techniques for the total assessment of agricultural 
targets. The LACIE procedure includes the collec- 
tion of a limited amount of meteorological data and 
the deployment of ground-truth reflectance panels 
adjacent to the agricultural targets. These targets 
are overflown by a multitude of sensor systems in- 
cluding the 11-channel modular multispectral scan- 
ner, which generated the data used in this report. 
The data used were collected over an area in Hand 
County, South Dakota, situated at 451 meters 
(1480 feet) above sea level. Flight dates were Sep- 
tember 25, 1975; July 8, 1976; May 15, 1978; and 
June 2, 1978. The data for these dates have been 
designated as Groups I, II, III, and IV, respectively, 
where used in the tables. Five reflectance panels 
with nominal reflectances of 60, 28, 16, 8, and 5 per- 
cent (figs. 11 and 12), with dimensions of 6.1 by 
12.2 meters (20 by 40 feet), are deployed side by 
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FIGURE 9. — Nomograph of scanner parameters for clear con- 
dition; visibility = 15 kilometers. 



FIGURE 11. — ^Reflectance as a function of wavelength for 5-, 
8-, 16-, and 28-percent-reflectance panels (September 25, 
1975). 


I 1 1 L 1 J 

* .3 .4 .5 .6 .7 .8 .9 1.0 

Wavelength, /xm 

FIGURE 10. — Nomograph of scanner parameters for excep- 
tionally clear condition; visibility = 100 kilometers. 



FIGURE 12. — Reflectance as a function of center wavelength 
for 60-percent-reflectance panel (September 25, 1975). 
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CaSft PROGRRM 


side. The data collection altitude is 457.2 to 609.6 
meters (1500 to 2000 feet), so that, with the 2.5- 
milliradian instantaneous field of view (IFOV) of 
the M^S sensor, between 8 and 10 data points 
(resolution elements or “reselms”) are obtained 
each scan line for 4 or 5 scans as the scanner over- 
flies a panel. Each reselm is digitized into 8 bits and 
recorded on tape so that the panel data are located 
in the printout through inspection by finding 
clusters of similar numbers in a 4 by 8 or a 5 by 10 
matrix. (See data in fig. 13.) These matrices are cor- 
related with specific panels by virtue of the corre- 
spondence of higher numbers with higher reflec- 
tances. When this correlation of signals with reflec- 
tance panels has been accomplished, the next step 
in computing NE Ap is to obtain a noise figure or 
signal. 

Noise is measured by noting the variations in 
the system output as the scanner is subjected to a 
uniform target. A convenient target of 8-reselm size 
is provided by the internal calibration or reference 
source, which is viewed once each scan. The noise 
is considered to be the standard deviation for the 
scanner output when the calibration source is 
viewed. Sufficient information is then available, 
using equation (23), to calculate NE Ap from direct 
measurements. 



where P 2 and p^ are the reflectances of any two 
panels, V 2 and are the corresponding average 
counts, and is the standard deviation of the 
counts while the scanner is viewing the calibration 
source. 

Each of the 10 visible channels for the M^S was 
analyzed by plotting the value of the reflectance of 
the panel taken at the center of the half-power 
bandwidths, as determined from laboratory cali- 
bration reports.2»3 By taking the low blackbody 
reference source as “zero” reflectance, a straight 
line was obtained as shown in figure 14. If the 
reflectance change and signal change between each 
panel are used, the result is a second straight line 
the slope of which is the same as the first but which 

2R. W. Nuss and J. R. Woodfill, Project ERAP Modular 
Multispectral Scanner Data Evaluation Report, JSC-09257, Oct. 
1974. 

^D. Hawley, M^S Engineering Baseline and Calibration Test 
Report, NASA-LEC 4969, 1975. 
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FIGURE 13. — Sample printout of M^S data showing clustering 
method used for locating reflectance panel data. 
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AP, percent 

FIGURE 14. — ^Plot of scene signal voltage as a function of rela- 
tive reflectance for channel 7 (September 25, 1975). 
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also passes through the zero signal-to-noise point 
(fig. 15). The slope of the reflectance as a function 
of the signal-to-noise ratio is the NE Ap for the 
channel. 

It would be profitable at this point to discuss 
why the line passes through zero in one case but 
not the other. First, it can be seen from the graph in 
figure IS that the line not passing through zero can 
be corrected to zero reflectance by a simple subtrac- 
tion of a constant offset count. The source of the 
offset term can be understood by reference to equa- 
tion (17), which, upon integration over wavelength, 
provides scene radiance. 

+ < 40 ) 

In equation (40), note that, for p = 0 percent, the 
scene radiance observed by the scanner is equal to 
the backscatter term Ng. Letting the responsivity of 
the scanner system be /? = 1^ /N^, then, in terms of 
counts. 


A schematic of the signal as a function of time is 
shown in figure 16, in which represents the zero- 
offset term. Thus, a panel-to-panel measurement is 
simpler than a panel-to-zero-radiance measurement 
because the backscatter term disappears. 


Comparison of Measured 
and Calculated NE Ap 

Values of NE Ap for the data acquired for the 
four flight dates of September 1975, July 1976, May 
1978, and June 1978 are recorded in table IV. In- 
cluded in this table are the Perry Moon values cal- 
culated from equation (28), the Angstrdm values 
calculated using equation (33), and the values ex- 
tracted using the nomograph discussed previously. 
In the P. Moon values, scene irradiance H is taken 
from table I, and and S/N are sensor perform- 
ance parameters taken from laboratory tests. As 
previously stated, the P. Moon values are calculated 
assuming an air mass of 2.0, which corresponds to a 
Sun angle of 60° measured from the zenith. The P. 
Moon value, then, is a rough calculation meant to 
produce a sensitivity number with regard to possi- 
ble, but not necessarily actual, conditions. In the 
Angstrdm values, NER is taken from laboratory 
data, and Hj, comes from the use of the Thekaekara 



FIGURE 15. — ^Plot of relative reflectance as a function of 
scene signal voltage for channel 3 (September 25, 1975). 



Vgg = Low BB count 

Vs * Scene signal 
(uncorrected) 
Vs' = Scene signal 
(corrected) 

Vg * 0-percent offset 


FIGURE 16. — Scanner video trace. 


table for the extraterrestrial solar irradiance (table 
V). 

An inspection of table IV reveals that for data 
Groups I and III, the P. Moon data more closely ap- 
proximate the measured values than do the 
Angstrom data. The Angstrom values are closest 
for data Group IV, whereas the nomograph pro- 
vides the best answers for data Groups I and II. 
Although the nomograph answers are not the best 
for Groups III and IV, there is more leeway for 
altering the numbers through the collection and use 
of more precise meteorological data. This adjust- 
ment would bring the measured and calculated 
numbers into better agreement. 

One large source of error in calculating NE Ap 
for the M^S is brought about by the separate optical 
paths of the ground scene and the calibration lamp. 
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II 


The result is that the NER term used in the calcula- 
tion is strongly influenced by the condition of the 
calibration optics, whereas the measured value of 
NE Ap is not influenced at all by these optics. This 
result indicates a need to monitor the condition of 
these optics periodically and suggests that NER 
measurements would be better made using external 
sources. 

Since making real-time measurements of NE Ap 
in an operational sense is generally impractical, the 
question arises as to the relative merits of the 
various methods of calculating NE Ap. The answer 
depends on several considerations including the 
use that will be made of the data. Obviously, the P. 
Moon values are the correct order of magnitude 
and are easily accomplished, but the other calcula- 
tions are capable of greater precision. 

Generally, sensitivity data are used for “ball- 
park” estimates of performance and it appears the 
P. Moon values accomplish this satisfactorily. 
However, it would be constructive at this point to 
consider the magnitude of variations in NE Ap that 
might occur because of such factors as time of year, 
time of day, latitude, and atmospheric conditions. 
This estimate can best be accomplished with the 
aid of the nomographs discussed previously. 

Figure 7 indicates that at or near the summer 
solstice (June 21), the air mass m is nearly equal to 
1.0 whatever the latitude, and this holds true for 
almost 2 hours before or after noon. Proceeding 
with this calculation for the data of June 2, 1978, 
channel 4 of the M^S, assuming clear conditions, 
one obtains for NE Ap a value of about 0.16 per- 
cent, which compares with a measured value of 
0.21 percent. The alteration of this number for hazy 
or exceptionally clear conditions is insignificant, 
but only near June. The situation can be vastly 
different, however, at the winter solstice (Decem- 
ber 22). Again choosing channel 4, it can be seen 
that air mass m may vary from about 1.7 at 30° 
latitude (at noon) to about 2.75 at 45° latitude. This 
variance results in a variation in the nomograph 
number from 1.6 to 5.3. These numbers in turn 
translate to a change in NE Ap of 0.31 to 1.0 per- 
cent. Additional changes will be brought about if 
the time is different from noon (as much as 42 per- 
cent in 2 hours), and this for a clear day. If haze is 
considered, the values for NE Ap can range from 
0.49 percent to nearly 2.0 percent at the two 
latitudes. For clarification, these values are given in 
table VIII. 


Table VIIL — Possible Variations in Calculated 
Values of NE lip for Channel 4 of the M^S 


NEip, 

percent 

Local time 

Month 

Latitude, deg 

Atmospheric 

condition 

0.15 

12m.a 

June 

30 

H‘aze to clear 

.16 

12 m. 

June 

45 

Haze to clear 

.31 

12 m. 

Dec. 

30 

Clear 

.49 

2 p.m. 

Dec. 

30 

Clear 

1,0 

12 m. 

Dec. 

45 

Clear 

2.0 

12 m. 

Dec. 

45 

Haze 

.95 

2 p.m. 

Dec. 

30 

Haze 

1.2 

12 m. 

Oct. 

45 

Haze 

.61 

12 m. 

Oct. 

45 

Clear 


®12 m. — noon. 


It should be argued that in agricultural surveys, 
data are rarely taken beyond the month of October. 
Then, if a limitation is also placed on atmospheric 
conditions so that data are acquired under clear air 
conditions, a figure of 0.61 percent is found for 
noon at 45° latitude. The corresponding P. Moon 
value is 0.39, which is 36 percent less. 


Conclusion 

The purpose of calculating noise equivalent 
reflectance is generally to provide an approximate 
value to a potential data user for use in estimating 
sensor performance under field conditions. Con- 
sidering the wide range of possible variations in 
noise equivalent reflectance that can be attributable 
to known conditions, as well as the influence of 
unknown, unpredictable, or time-varying condi- 
tions, it is believed that the P. Moon calculation 
provides the most convenient and usable indicator 
of performance. One outstanding reason for this 
conclusion is the simplicity in defining a standard 
procedure as opposed to defining standard condi- 
tions, with artificial compromises or constraining 
factors that would significantly complicate the pro- 
cedure. 

Situations may arise wherein the success of a 
particular experiment may be marginally depend- 
ent on the sensitivity of a scanner system under 
certain conditions. In such situations, a direct 
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measurement using reflectance panels is recom- 
mended. Otherwise, the nomograph described in 
this report should be used as the best indicator of 
system sensitivity. 


Lyndon B. Johnson Space Center 
National Aeronautics and Space Administration 
Houston, Texas, July 16, 1979 
953-36-00-00-72 
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Appendix A 

Meteorological Effects on Radiation 


An airborne scanner system collects radiant 
energy that has been reflected or emitted by a target 
on the surface of the Earth. The attenuation of 
solar radiation through the atmosphere from its ex- 
traterrestrial value is determined by the state of the 
atmosphere at any given time. The physical condi- 
tions of the airmass can be highly variable depend- 
ing on temperature, humidity, dust content, etc. An 
exact model for radiation attenuation has not been 
developed, but approximations can be made in 
which the gross, average irradiance at the surface 
can be estimated. An extensive study on measured 
atmospheric quantities and radiation transfer 
models has been made (refs. 5, 7, 8, and 9). Only 
those features of interest to the current report will 
be discussed here, and the discussion is not to be 
considered as an exhaustive account. 

In figure 17, consider a beam of parallel direct 
solar extraterrestrial spectral irradiance (watts 
per square centimeter per micrometer) at the top of 
a stratified atmosphere. For a parallel beam trans- 
mitted a distance dz, the irradiance will decrease by 
and assuming no radiant emittance within the 
distance dz, Schwarzchild’s equation gives 


H. 


-k^dz 


(Al) 


where is the value of the irradiance at the top of 
the layer and is the volume extinction coeffi- 
cient (attenuation coefficient including both ab- 
sorption and scattering effects) (per kilometer). In- 
tegrating equation (Al) gives 




-fo k^dz 


(A2) 


and for a surface FA at some angle to the zenith, 


\ cos c 

However, dx = sec 9 dz, 


-/o *x 


dx 




(A3) 


(A4) 


where 


-I 


(sec6)z 


k^dz 



FIGURE 17. — Geometric relationships of solar irradiance 
quantities used in Schwarzchild’s equation. 


equals optical depth or thickness of the atmo- 
sphere. The quantity sec 0 is normally called the air 
mass m of the atmosphere. This can be seen by con- 
verting equation (A3) in terms of a mass extinction 
coefficient where p equals mass den- 

sity. Letting dm equal mass element per unit area 
along a column of air in direction dx, one can write 


dm = p dx 

- p sec 0 dz 


(A5) 


Integrating and letting equal the mass per unit 
area in a vertical column. 


m - (sec0)m, 


0 


(A6) 


Hence, sec 9 = m/m^ = m^. Equation (A4) can 
now be written as 


^x = ^Xq cos 0Q e 


(A7) 
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In general, the optical depth is determined by 
both absorption and scattering. Angstrom’s at- 
mospheric model for the visible spectrum (ref. 7) 
considers an average atmospheric state such that 


r 




+ T + T 


(A8) 


where r^is due to Rayleigh (molecular) scatter- 
ing including corrections for refraction and 
anisotropy of the air molecules, haze, or dust scat- 
tering, and T,„is the absorption due to water vapor, 
ozone, carbon dioxide, and other gases. Sometimes, 
the quantity is written as e~'”r^ — 

where = e”'’" and is termed the transmissivity of 
the atmosphere. For analysis of NE Ap, Ang- 
strom’s model is used such that 




(a9) 


the Earth’s surface, and that is the irradiance due to 
multiple scattering of light within the atmosphere. 
As might be expected, this diffuse sky irradiance 
is highly variable and no satisfactory atmospheric 
model fits experimental data taken at various loca- 
tions. Models of 5^ have been determined using 
scattering functions in the radiation transfer equa- 
tion (ref. 5), and the concept is expressed by the 
simple formula 

— = kco%Qf. (AlO) 

^x 

where k = k(Tjy,Tg,m) and is the direct solar ir- 
radiance at the surface. Hence, the total irradiance 
of a target scene is given by 


= ^x. e 




+ 5 , 


(All) 


where (Rayleigh extinction coefficient) 

= 0.00889X~* ®® for X in micrometers 
Tg = Gg (dust, haze extinction coeffi- 
cient) 

Gg = 2.303R(2X)”“ for X in micrometers 
= G^ (absorption extinction coeffi- 
cient) 

To fit the modular multispectral scanner data for 
latitude 45° N (the position of the panels), the 
average turbidity B for this latitude and absorption 
are foimd in charts given in references 5 and 8. 
Finally, for Gg, the normal value of a for haze scat- 
tering is taken as 1.5. As the haze particles change 
in size and number per unit volume, a will lie be- 
tween 0 a < 4, where for 1 -micrometer size par- 
ticles, a ®= 0.5. 

One other atmospheric effect must be included 
to obtain the total spectral irradiance on a target on 


Another useful atmospheric model used in find- 
ing optical depths in terms of the horizontal visual 
range is found in reference 10. In this model, the 
aerosol (haze) optical depth Tg is related through 
the turbidity of the atmosphere and the volume ex- 
tinction coefficient to the horizontal visual range. 
The calculated Rayleigh scattering and aerosol opti- 
cal depths as a function of wavelength and visual 
range are shown in figures 5 and 6. These were used 
in calculating the range of values for the ratio of 
system parameters [//g/n-][NE Ap/NER] shown in 
figures 8, 9, and 10. In addition, the ratio of sky irra- 
diance to direct solar irradiance was taken as 




^X 


cos 0Q 

(A12) 

_ 0^ 
m 

(A13) 
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Appendix B — Symbols 


A, 

collecting area, cm^ 

Hsx 

total irradiance of target scene 


target area, cm^ 

Ht 

total target irradiance, W-cm”^ 


aerosol attenuation coefficient 

Hax 

solar irradiance, W«m“^»/Ltm~^ 


Rayleigh extinction coefficient 


(table I) 


water-vapor absorption coefficient 

Hx 

spectral irradiance, W«cm“^* 
^im“^ 

B 

turbidity coefficient, 0.05 ^ B ^ 



0.2 


spectral normal irradiance 

Z)(0-X) 

percentage of solar constant asso- 
ciated with wavelengths shorter 

^0 

extraterrestrial solar irradiance; 
solar constant, W-m“^ (table I) 


than \ (table V) 

h 

hemisphere 

dA 

element of area, cm^ 

qm-- 0.68 


dfix 

decrease in irradiance 


volume extinction coefficient. 

dm 

mass element per unit area in 


km“^ 


direction dx 

m 

absolute air mass, unitless; ap- 

dPi 

power incident on surface 


proximately equal to secant 9 for 

dPr 

power reflected from surface 


low altitudes and for 0° ^ ^ 60° 

dx = sec 0 dz 


= m/njQ 

relative air mass, unitless 

dz 

beam transmission distance 

= sec^ 


d€l 

projected solid angle 


mass per unit area in vertical col- 

E{\) 

solar spectral irradiance averaged 

N 

umn 


over small bandwidth centered at 
\, W.m“^»/xm“‘ (table V) 

radiance (power radiated), 
W-cm“^* sr“* 

£(0-X) 

area under solar spectral irra- 
diance curve in wavelength range 

Ha 

backscatter radiance, W«cm“^- 
sr“* 


0to\,W.m-2 (table V) 

Hci\ 

radiance produced by calibration 

H 

irradiance (power density), 


lamp 


W*cm“^ 

Hp 

target radiance, computed from ta- 

ffa 

irradiance, AngstrSm’s model, 


ble! (P. Moon), W*cm”^*sr~^ 


W-cm-2 

Hr 

path radiance due to scattering 

Hb 

total irradiance from all sources, 

W.cm”2 


into beam from surrounding 
targets 

He 

extraterrestrial solar irradiance. 
Turner’s model 

HrX 

path spectral radiance due to scat- 
tering of radiation into the beam 

Heb 

extraterrestrial solar irradiance in 


from surrounding ground targets 

given spectral band, W -cm”^ 

Hs 

scene radiance, W«cm“^*sr~* 

^ex 

extraterrestrial spectral solar irra- 

Hr 

target radiance, W«cm“^*sr“* 

diance, W-cm“^»/Ltm”^ 

Hrx 

scene or target spectral radiance 

Hex 

parallel direct solar extraterrestrial 
spectral irradiance, W»cm“^« 


spectral radiance, W-cm'^-sr™'- 
/xm“* 



\(cnl) 

calibration source spectral ra- 

H, 

ground irradiance, W -cm”^ 

diance 

H,x 

ground spectral irradiance 

NER 

noise equivalent radiance. 

Hp 

target irradiance from table I (P. 


W-cm”^-sr“^ 


Moon), W-cm*“^ 

NER, 

spectral noise equivalent radiance 
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NE Ap 

Pi 

Pm 

Pr 

m 

q = ^ 

R = V,/N, 

P. 

S 

Sb 

S, 

S, 

S/N = VJV„ 
T 
V 
Va 

^BB 

»;al 

Vn 

Vs 

v; 

Vl,V2 

x.y.z 

“px= ^\/p 
A 

^N, 

AiV, 

Ap 

e 

X 

^ c 
>^0 


noise equivalent reflectance, unit- 
less 

power incident, W 

power input, W 

power reflected, W 

surface at some angle to the zenith 

(fig. 17) 

transmissivity of atmosphere 
scanner system responsivity, 
V*(W.cm“^*sr~*)”^ (eq. (41)) 
visual range 

sky irradiance, W»cm“^ 
diffuse sky irradiance in given 
spectral band pass 
sky irradiance for Turner’s equa- 
tion, W*cm"^ 

diffuse or sky spectral irradiance 
for Angstrom’s equation, 
W*cm”^* pm“‘ 
signal-to-noise ratio, unitless 
time, hours 
emf, V 

zero-percent offset emf, V 
emf, low blackbody count, V 
calibration emf, V 
noise emf, V 

scene signal emf, uncorrected, V 
scene signal emf, corrected, V 
scanner output emf from targets 1 
and2,V 

rectangular Cartesian coordinates 
mass extinction coefficient 
change 

change in scene radiance 
change in spectral radiance 
change in scene reflectance 
nadir angle, deg 

nadir angle of reflected solar en- 
ergy, deg 

nadir angle of incident solar en- 
ergy, deg 
wavelength, pm 
center wavelength, pm 
a selected wavelength, pm 


P 

p' 

Peal 

Peq 

Pn 

Ps 

PQ 


~ ^R 


<i>r 

a 


Subscripts: 

a 

BB 

b 

c 

cal 

e 

eb 

eq 

ex 

g 

h 

i 

in 

n 

P 

R 

r 

s 

t 

V 

w 

\ 


total reflectance; mass density (eq. 
(A5)) 

partial reflectance, sr"^ 
calibration reflectance 
equivalent scene reflectance 
noise reflectance 
signal reflectance 
a selected value of reflectance 
optical depth (optical thickness) 
of atmosphere 
aerosol (haze) optical depth 
total Rayleigh optical depth 
absorption due to water vapor and 
other gases (eq. (A8)) 

azimuth angle of reflected solar 
energy, deg 

solid angle of reflectance, sr 
solid viewing angle, instantaneous 
field of view,sr 


aerosol 

blackbody 

band (band pass) 

collecting (collector); center 

calibration 

extraterrestrial 

extraterrestrial in given band 

equivalent 

extraterrestrial 

ground 

hemisphere 

incident 

input 

noise 

P. Moon 

Rayleigh 

reflected; relative 
scene (signal) 
target; Turner 
visual 

water vapor; water vapor and 

other gases 

wavelength 
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p 

0 


mass density 
reference (selected) 


Superscripts: 

m 

a 


relative air mass 

Angstrom’s exponent, 0 < a ^ 4 
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